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xi?£um  gas  temperature  calculations  ha^H&eetr-crJfnpleted 
for  the  7/8  inch  I.D.  test  section.  The  results  of  these 
calculations  are  showru-jn  Fig.  1 . 

2.)*i*Coinbustion  tests  were  completed  for  Cetane  and  Mil  2190 
/TEP  in  the  7/8  inch  l.D.  test  section.  Two  methods  of 
fuel  introduction  were  used,  microfog  and  a  wick  saturated 
with  the  test  fluid.  The  results  of  these  tests  are 
summarized^in  Figs.  2  and  3. 

By  employing  SIT  data  for  Mil  2190  as  given  in  Ref.  4  a 
combustion  limit  curve  has  been  calculated  for  this  fluid 
in  the  7/8  inch  I.D.  test  section.  A  comparison  of  the 
calculated  curve  with  experimental  results  is  shown  in 
Fig.  4.  The  calculated  curve  is  somewhat  more  conservative 
than  the  experimental  results  as  would  be  expected. 

4.>  The  combustion  test  modifications  completed  for 

the  3/8  inch  and  1  1/2  inch  I.D.  test  sections . 

-5-r^  An  analysis  for  the  case  where  the  charging  process  is.e<e^ 
governed  only  by  entrance  and  throttling  losses  has  been7 
made  for  the  liquid  charging  case.  Expressions  for  predict¬ 
ing  the  pressure  vs .  time  curve  and  the  maximum  pressure 
rise  rate  were  derived  and  show  good  agreement  with  test 
results  fiD V/7S9 IJb  Hn&09a»d!4D  HHtfo 
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II.  Maximum  Gas  Temperature  Calculations 

For  the  maximum  gas  temperature  calculations  the  analysis 
given  in  Ref.  1  was  employed.  The  convection  correlation  used  in 
the  analysis  was  that  derived  in  Ref.  2  for  combined  free  and  forced 
convection.  For  these  calculations  the  following  supplementary 
assumptions  were  made: 

1. )  The  effect  of  axial  heat  conduction  is  negligable. 

2. )  The  maximum  gas  temperature  is  attained  at  the  center  of  the 

fluid  increment  extending  one  foot  upstream  of  the  end  plate 
prior  to  the  compression  process. 

As  pointed  out  previously  (Refe.  1  and  2)  both  these  assumptions 
are  consistent  with  the  experimental  results  obtained  to  date. 

The  calculations  assume  an  initial  gas  and  pipe  wall  tempera¬ 
ture  of  80°f  with  the  initial  pressure  in  the  downstream  pipe  being 
14.2  psia  which  corresponds  to  the  average  conditions  in  this 
laboratory. 

The  calculation  procedure  was  as  follows: 

1. )  A  pressure  rise  rate  was  selected. 

2. )  The  calculation  proceeded  step  by  step  with  the  gas  temperature 

at  the  center  of  the  increment  being  recorded  as  a  function  of 
gas  pressure. 

3. )  This  procedure  was  continued  until  a  point  was  reached  where 

the  gas  temperature  began  to  decrease  as  the  pressure  increased. 
The  maximum  temperature  attained  at  any  pressure  higher  than 
the  pressure  at  which  the  gas  temperature  began  to  decline  was 
taken  as  being  the  maximum  temperature  given  by  calculation. 

By  continuing  this  process  for  a  variety  of  pressure  rise  rates 
eventually  enough  data  were  accumulated  to  allow  a  plot  of  constant 
temperature  lines  with  respect  to  pressure  rise  and  maximum^ pressure. 
The  results  are  shown  in  Fig.  1  for  the  gas  temperature  range  between 
250  and  900°F,  which  is  the  most  critical  temperature  region  in 
regard  to  spontaneous  ignition. 

The  interpolation  for  plotting  constant  lines  was  alone  on  semi¬ 
log  paper  in  order  to  take  into  account  the  exponential  change  in 
temperature  with  a  linear  change  in  pressure. 
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II.  Combustion  Tests 

2.1  Test  Procedure 

The  7/8  inch  I.D.  test  section  was  employed  for  these  tests, 
the  test  section  being  5  feet  lone  and  in  a  horizontal  position. 

The  procedure  for  conducting  the  miscofog  tests  has  already 
been  outlined  in  Ref. 3. 

For  the  wick  tests  a  small  piece  of  fiberglass  cloth  was 
attached  to  a  post  along  the  centerline  of  the  test  section,  extend¬ 
ing  from  the  end  plate  to  1/2  inch  upstream.  Prior  to  a  test  this 
wick  was  saturated  with  a  measured  quantity  of  the  test  fluid, 
applied  with  a  hypodermic  needle.  Whether  combustion  occurred  or 
not  the  wick  was  replaced  after  each  test.  Ignition  was  detected 
both  by  the  photocell  and  by  visual  observation  of  the  wick  follow¬ 
ing  a  test.  In  all  instances  these  two  methods  of  ignition  were 
in  complete  agreement. 

For  both  forms  of  fluid  introduction  a  test  condition  was  only 
noted  as  a  point  of  no  ignition  if  ignition  failed  to  occur  for  4 
tests. 

2.2  Test  Results 

The  results  of  the  Cetane  and  Mil  2190  TEP  tests  are  plotted 
in  Figs.  2  and  3  respectively.  For  bothof  these  figures  microfog 
ignition  is  denoted  by  the  symbol  (#)  and  ignition  with  the  wick 
by  (A)  ,  points  of  no  ignition  are  indicated  by  the  same  symbols 
( 0,A )  not  blacked  in. 

The  line  issuing  from  the  origin  of  the  coordinate  system 
and  slanting  upward  and  to  the  right  is  the  pressure  rise  rate 
limit  line.  The  line  denotes  the  maximum  attainable  pressure  rise 
rate  for  a  given  final  pressure  with  the  present  test  configuration. 

The  microfog  tests  with  Cetane  were  conducted  with  an  overall 
air-fuel  ratio  of  59.  Further  tests  were  made  using  higher  air -fuel 
ratios  with  no  change  in  the  results  shown  in  Fig.  2  until  the  air 
fuel  ratio  reached  180  above  which  slightly  higher  final  pressures 
and  pressure  rise  rates  were  necessary  for  ignition.  The  microfog 
tests  with  Mil  2190  TEP  were  conducted  with  an  overall  air-fuel 
ratio  of  113  which  was  the  lowest  that  could  be  attained  for  this 
fluid  with  the  present  microfog  unit.  In  this  case  detectable 
variations  in  the  combustion  limit  curve  were  noted  when  the  air 
fuel  ratio  reached  200. 
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2.3  Discussion 

The  fact  that  the  mode  of  fluid  introduction  has  a  different 
effect  on  different  pressure  ranges  points  out  the  necessity  of 
including  both  microfog  and  wick  tests  in  order  to  find  the  most 
conservative  results.  Whether  some  other  form  of  fluid  introduction 
would  give  more  conservative  results  remains  an  open  question. 

The  physical  reason  for  such  behavior  is  not  yet  apparent  and  it  is 
hoped  that  the  work  to  be  undertaken  in  the  next  few  months  will 
help  to  clarify  this  picture  somewhat. 

Although  the  results  of  the  present  tests  with  Mil  2190  TEP  and 
those  completed  by  Electric  Boat  Division  as  replotted  in  Ref.  4 
are  not  strictly  comparable  due  to  difference  in  pipe  diameters  and 
maximum  pressure  it  is  of  interest  to  note  that  the  minimum  pressure 
rise  rates  for  combustion  are  of  about  the  same  order  of  magnitude. 
The  present  results  are  not  as  conservative  which  would  be  expected 
due  to  the  higher  heat  losses  in  the  smaller  pipe  and  the  increase 
in  combustion  resistance  at  lower  pressures  as  evidenced  by  the 
Self  Ignition  Temperature  (S.I.T.)  curve  for  Mil  2190  TEP. 

Employing  the  data  on  the  S.I.T.  of  Mil  2190  TEP  at  various 
pressures  as  given  in  Ref.  4,  a  theoretical  combustion  limit  curve 
was  determined  for  this  fluid.  For  this  determination  the  isotherms 
of  Fig.  1  were  employed,  assuming  that  combustion  would  occur  at  a 
given  pressure  when  the  gas  temperature  reached  the  S.I.T.  This 
result  is  shown  in  Fig.  4  where  the  experimentally  determined 
combustion  limit  curve  for  the  wick  tests  is  also  included  for  com¬ 
parison. 

As  would  be  expected,  the  theoretical  curve  is  somewhat  more 
conservative  than  the  experimental  one.  This  is  due  to  the  fact 
that  the  ignition  lag  is  entirely  neglected  in  the  analysis  and 
also  because  the  calculated  temperature  given  by  the  analysis  of 
Ref.  1  have  always  been  above  the  measured  temperature. 

III.  Liquid  Changing  Study 

3.1  Analysis 


An  analysis  has  been  developed  and  is  presented  in  the  Appendix, 
for  the  case  where  the  charging  process  is  governed  only  by  entrance 
and  throttling  losses.  Two  cases  were  considered:  (1)  the  air 
column  is  compressed  isentropically  (2)  the  air  is  compressed  iso¬ 
thermal  ly. 


The  assumptions  made  throughout  the  analysis  are  as  follows: 
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(A)  Air  Column 

(1)  A  linearized  representation  of  the  isentropic 
compression  process  is  used. 

(2)  The  ratio  of  specific  heats,  K,  is  assumed 
constant  and  equal  to  1.40. 

(B)  Charging  Liquid 

(1)  The  density  of  the  liquid  remains  constant. 

(2)  The  pressure  and  velocity  of  the  liquid  are 
assumed  uniform  across  any  cross  section  except 
in  the  immediate  vicinity  of  the  orifice  and  the 
entrance  to  the  upstream  pipe. 

(3)  The  effects  of  fluid  acceleration  are  neglected. 

(4)  Friction  forces  exerted  on  the  liquid  by  the  pipe 
walls  are  neglected. 

(5)  The  coefficient  of  discharge  of  the  orifice  is 
assumed  equal  to  1.0. 

(6)  Quasi-steady  state  conditions  are  assumed  at  the 
orifice. 

(7)  Changes  in  elevation  are  neglected. 

(8)  The  valve  opens  instantaneously. 

(C)  Air  in  the  Storage  Tank 

(1)  The  air  is  expanded  isentropically . 

(2)  During  the  process  the  change  in  air  volume  in  the 
tank  is  negligible. 

The  analysis  provides  a  method  for  predicting  the  pressure  vs. 
tiihe  curve  and  the  maximum  pressure  rise  rate. 

The  equations  for  predicting  the  pressure  vs.  time  curve  are: 

(in  the  notation  given  in  the  Appendix) 
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where : 


Ci  ~  (l  -  £,)  +  K  (\ -  §0)fa 
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The  maximum  dimensionless  pressure  rise  rate  is  given  by: 

=  0.287^  =  «287(!-S.) 
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a  p4> 

1 
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(3) 


For  the  linearized  isentropic  case  with  gas  pressures  in  the  range 
15  to  2000  psia,  K  ■  1.40,  ■  .95.  For  the  isothermal  case, 

K  -  1, 


It  is  seen  from  Gq.  3  that  the  maximum  dimensionless 
rise  rate  is  a  function  of  the  dimensionless  number, 

y 

O  L  I  ?  Wn  I  12. 

N  =  (I  -  f )  te 


pressure 


pj.  X, 


Z  Pc. 


W-H) 


which  depends  upon: 

(1)  the  driving  pressure,  P^. 

(2)  the  downstream  pipe  length,  X2- 

ftp 

(3)  the  pipe  area  to  orifice  area  ratio,  jp 

(4)  the  nature  of' the  compression  process,  whether  it  be 
isentropic  or  isothermal. 


A  comparison  between  the  theoretical  J  t  vs .  time  curve 
and  the  experimental  pressure  trace  is  shown  in  Fig.  5. 

In  Fig.  6  the  maximum  normalized  pressure  rise  rate, 
is  plotted  versus  the  number,  LL  ,  with  the  reference  time 
unit  To  equal  to  1  sec.  *■<> 

3.2  Tests 


A  series  of  tests  were  run  in  order  to  verify  the  analysis. 
The  following  parameters  were  varied: 

(1)  Driving  pressure  -  teats  were  run  at  600,  1200,  and 
2000  psig . 

(2)  Orifice  size  -  diameters  of  .070,  .052,  and  .035  in. 
were  used. 

(3)  The  length  of  the  downstream  pipe  section  lengths  of 
60  and  33  in.  were  used. 

For  all  tests  the  following  parameters  were  held  constant: 

(1)  Upstream  pipe  length  ■  184  in. 

(2)  Valve  opening  time  -  16  m  sec. 

(3)  Pipe  inside  diameter  -  .375  in. 

3.3  Results  and  Discussion 


Reasonable  correlation  between  theoretical  and  experimental 
results  was  obtained  as  shown  in  Figs.  5  and  6.  A  comparison  of 
the  experimental  and  theoretical  results  shows  that  the  compression 
of  the  air  column  is  somewhere  between  isentropic  and  isothermal. 

The  parameter  N  in  Eq.  3  is  of  interest  since  it  reveals  the 
variables  influencing  the  maximum  pressure  rise  rate.  It  is  seen 
that  the  maximum  pressure  rise  rate  is  directly  proportional  to  the 
driving  pressure  raised  to  the  3/2  power,  inversely  proportion  to 
the  downstream  pipe  length,  and  varies  in  an  inverse  manner  with  the 
pipe  Co ‘.orifice  area  ratio.  It  is  also  influenced  by  the  nature  of 
the  compression  process,  whether  it  be  isentropic  or  Isothermal. 

It  is  important  to  note  that  the  analysis  developed  assumes  that 
the  charging  process  is  governed  only  by  entrances  and  throttling 
losses.  It  remains  to  develop  a  method  for  determining  under  what 
conditions  these  assumptions  may  be  used. 
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Future  Work 


1.)  Combustion  tests  will  be  continued  using  Mil  2190  TEP 
in  the  3/8  inch  I.D.  test  section,  both  microfog  and 
wick  tests  will  be  included. 


2. )  A  series  of  napthenic  neutral  oils  have  been  obtained 

with  flash  points  varying  between  300  and  500<>f.  These 
oils  will  be  tested  to  determine  the  relative  advantage, 
if  any,  of  employing  fluids  with  a  higher  flash  point. 

3. )  Maximum  temperature  calculations  will  be  made  for  the 

3/8  inch  I.D.  test  section  in  order  to  afford  a  comparison 
with  the  measured  results  similar  to  that  given  in  this 
report. 

4. )  With  regard  to  the  liquid  charging  tests  an  effort  will  be 

made  to  ascertain  under  what  condition  the  included 
analysis  may  be  used. 

5. )  A  study  will  be  made  to  determine  the  maximum  pressures 

attainable  if  combustion  occurs  during  the  liquid  charging 
process . 
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Control  Volumes  Employed  in  Analysis 


Appendix 

Analysis  of  Liquid  Charging  Process 

A-l  The  analysis  for  the  case  where  the  process  is  governed  only 
by  entrance  and  throttling  losses  is  presented  with  reference 
to  Fig .  7 . 

A-2  Notation 

A0  =  cross  sectional  area  of  orifice,  in2. 

Ap  “  cross  sectional  area  of  pipe,  in2. 

At  m  cross  sectional  area  of  storage  tank,  in2. 

cK  =  differential  area  vector,  in3. 

Cv  *  specific  heat  at  constant  volume  for  air.  Btu/lb. 

C] ,  C2,  C3  •••  =  dimensionless  constants  defined  throughout 
the  analysis  for  ease  of  manipulation. 

e  *  internal  energy,  in-lbf/lbm. 

Ff  ■  friction  force,  lb. 

p  =  external  force  acting  on  a  control  volume,  lb. 

2C  =  gravitational  constant,  lbm  in 
CF  TFT  sec2. 

g  =  acceleration  of  gravity,  in/sec2. 

h  =  enthalpy  in-lbf 
lbm 

k  ■  ratio  of  specific  heats  of  air. 

Pg0  “initial  pressure  of  air  in  downstream  pipe  section,  psia. 
Pg  -  pressure  of  the  air  in  the  downstream  pipe  section,  psia. 
PjJq  =  initial  air  pressure  in  the  storage  tank,  psia. 

Pj  •  air  pressure  in  storage  tank,  psia. 

pl»  p2»  p3»  • • •  “  pressure  acting  at  control  surfact?  l,  2f  3, 

...  respectively.  See  Fig.  7. 


li 


normalized  pressure  rise  rate,  sec.”l 

dimensionless  pressure  rise  rate. 

P shear  “  the  rate  at  which  shear  forces  dissipate  energy. 

Q  “  the  rate  at  which  heat  is  transferred  into  a  control  volume, 
t  =  time ,  sec . 

t  dimensionless  time  parameter,  *  £ 
tQ*  reference  time  unit,  sec. 
v  *  velocity  of  fluid  particles,  in/sec. 

V  =  volume,  in^. 

xl»  x2»  Yi,  Y2  are  defined  in  Fig.  7. 

x  *  system  coordinate  =  length  of  liquid  column  in  downstream  pipe 
section. 

y  *  distance  traveled  by  liquid  level  in  tank,  in. 
z  =  elevation  above  an  arbitrary  datum,  in. 

£  =  dimensionless  system  coordinate,  “  x 

X2 

point  about  which  the  function  (1  -  S)k  is  linearized. 

P  =  density  of  liquid,  lbm 

Tn3. 

density  of  gas,  lbm/in3. 
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A-3  Control  Volume  1-2 

The  momentum  equation  is  written  in  the  form 

2F  =  +  ffcv-clAi/ 

c.v.  c.  s. 


which  becomes 

-  r  Pz  A  P~  P-} A  p 


where  Fq  is  the  friction  force  exerted  on  the  liquid  by 
the  pipe  walls. 

Carrying  out  the  integration  and  differentiation  gives, 


Assuming  that  Fif  -  0  and  that  the  effects  of  fluid 
acceleration  are  negligible,  Eq.  2  becomes 


4  ^ 


3 
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A-4  Control  Volume  2-3 

The  momentum  equation  is  again  applied; 


Assuming  that  the  pressure  acts  uniformly  over  surface  3, 
the  above  reduces  to, 


+  ^3  Ap  '  %  Ap 


Neglecting  friction  and  assuming  quasi  -  steady  conditions 
gives , 


P3  -  Pa  = 


4 


For  incompressible  flow  the  continuity  relation  is 

4 1/1 '  4  a?  A° 1/3 

\4  =  Pf  d* 

A0  dt 


S’ 


V 


Substituting  Eq.  5  into  Eq.  4  gives, 


Pt.  ~  Pz  - 


4 


A-5  Control  Volume  3-4 

The  energy  equation  is  written  in  the  form, 

e  -  tL.r  *  jf/fl  e  dv  +  f[h  *  %](  M  7 

c,v5  C.  5, 


Assuming,  Q  -  P8hear  “  z  “  o,  the  internal  energy  of  the 
liquid  remains  constant  and  quasi-steady  state  conditions, 
Eq.  6  reduces  to, 


c.S. 


which  becomes 


vi 


The  continuity  relation  is 

/U  -  *  4$ 


Therefore  Eq.  8  reduces  to 


A -6  Control  Volume  4-5 

The  momentum  equation  is 


which  reduces  to 

P$  A?-  ft 


>D 


Neglecting  friction  and  assuming  the  effects  of  fluid 
acceleration  negligible  Eq.  10  becomes 


fc  -  ^4 


// 
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A -7  Control  Volume  5-6 

The  energy  equation  is  written  in  the  form 


/2 


The  following  assumptions  are  made: 

(1)  Kinetic  energy  of  the  gas  is  negligible. 

(2)  Internal  energy  of  the  liquid  remains  constant. 

(3) ’  Elevation  effects  are  negligible 


Q.  3  ** shear  =  ^ 


Eq.  12  then  reduces  to 


dA 


Carrying  out  the  integration  gives, 


P  (7.-  a ) 

. 


a 

ft  Mr 


t  3 


Assuming  an  ideal  gas  it  follows  that 

Uq  -  Cv 
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CjlP  =  JL 
K  /<-/ 


M 


The  gas  in  the  air  tank  is  assumed  to  expand  isentropically 
hence , 


Substituting  Eqs .  15  &  14  into  Eq.  13  and  differentiating 
gives; 


Neglecting  the  effects  of  fluid  acceleration  and  employing 
the  continuity  relation  A  qy  _  /)p  ?  Eq.  16 

becomes  t  dt 


ix 


Assuming  that  y  is  always  small  compared  to  Y,,  ie  Y,  +  q  ~  Y., 
Eq.  1 7  reduces  to,  ® 


& 


f  (d*  l* 


/f 


A-8  System  Differential  Equation 

Combining  Eqs .  18,  11,  9,  6,  and  3  gives 


The  air  column  in  the  downstream  pipe  section  is  assumed  to 
be  compressed  isentropically ,  hence 


X 


Eq.  20  is  substituted  into  Eq.  19  and  the  resulting  equation 
is  put  into  dimensionless  form  giving 


[U  )  IcH"/ 


Z./ 


where 


€  -  7C.X,  ±  ^  tot 


'to  -  c.  ■tiwi-e.  hvmX 


For  ease  of  manipulation,  and  C2  are  defined  as 

'"I  -Az^'T- 


and  Eq.  21  takes  the  form 


Z.Z 


The  initial  condition  for  the  system  is 
£  =  O  at  t*  -  O 
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A-9  Solution  of  the  Differential  Eq.  22 

The  function,  (1  is  linearized  by  the  use  of  a  Taylor 

Series  expansion,  neglecting  higher  order  terms 

$)*■ 

.f(£)~  +  ^  (£.!(£•  0+-  -  -  • 

■  Jr(£)~'  (|-0K-  <Ci-cf(S- f.) 

^5)  -  (i-gJWi-cfc  -  eU-cTf 


Eq.  24  is  then  of  the  form 

${%)  )K  ~  ~  C«  % 

0 -«„)**  <U-zfi 

C*  =  K  ( \  ~  ?„) 


Eq.  22  may  be  written  as 


l/E.  d? 


24 


Eq.  26  is  integrated  by  performing  the  following  transformation 


Let 


Ct>  -  Cq 


&■ 


Sj.c4,  Z^9tco0c\& 

-o,0 

d<f  ~  2si^&  ccz>Q 


Cq  Ci 


d?  ~  2  £i.  i£S8_  Jg, 

C< 


Eq.  26  now  becomes 


a  iff.  Si  de 

Co  S\v\b& 


-  sw\x© 


dt 


cX 


Z.\lc  Z.C, 


69 

Si 


c!  k 


2.1 


The  Boundary  conditions  of  Eq.  23  are: 


xiii 


From  Eq.  27  it  follows 

0/  a  Q  ojt  ^  ~  ® 

3 


0 .  t 


# 

t 


»  C> 


J.? 


Eq.  28  is  now  integrated  with  the  following  formula  (Ref.  5): 


r 


de  _ 


4  SI, 


1 


4.  cons'ta'AT 


Integrating  Eq.28  with  conditions  29  gives 

9 

.  &c, 


al/clc,  ( J4- 

Q 


c« 


e 


SlvA  t? 


i/5 


Putting  in  the  limits  and  employing  the  appropriate 
trigonometric  formulas  gives. 


C.ifZz. 


^4 


•  i/c. 


r-V'-ft  * 


<4 


& 


* 

.*£  m 


From  Eqs .  27,  25,  22,  21,  and  19,  it  follows  that 


V  c> 

/4«AA  O'  = 


^-c<f 


P»0 


<U/V\ 


]_£  L 

Eqs.  31  and  32  form  a  solution  from  which  the  Po0  ^ 
curve  may  be  constructed.  The  solution  assumes 
a  linearized  representation  of  the  function,  (1  -  ?)  , 
about  the  point  .  Depending  upon  the  range  of  pressure 
involved,  the  point  ^  will  vary.  A  value  of  -  0.95 
was  used  in  the  calculations  and  the  function  was 

represented  quite  well  in  the  range  o  <  C  =.  -3S" 
which  corresponds  to  a  pressure  range  for  the  air  column 
of  15  to  2000  psia . 


A-10  Isothermal  Compression  of  the  Air  Column 


The  solution  for  this  case  is  found  by  noting  that  for  an 
isothermal  compression  the  air  pressure  is  given  by: 


and  in  dimensionless  form 


(>-«) 


Comparing  Eq.  33  with  Eq.  20  it  follows  that  the  solution 
given  by  Eq.  31  and  32  reduces  to  the  isothermal  case  if 
K  -  1. 
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A-ll  Maximum  Pressure  Rise  Rate 

Recalling  Eq.  32  it  is  seen  that 

£}  -  &  ( repair j)  3  2. 

Poo  " 


■* 

Differentiating  with  respect  to  t  gives 


Z  si ^  B  coS  & 


c!& 

cit* 


d 2 
di* 


is  obtained  by  differentiating  Eq.  31  implicitly 


and  upon  substitution  into  Eq.  35  gives. 


0 /ff 


2  -su/i  &  col>6* 


Z  cos  & 


~\v\&  + 


2  4<1aa  & 

"Z. 


Cq 


■3  4* 


Noting  that 


«e‘f 


Z.'taA/'i  6  iM  £? 


-  «*ac.se  ^ 

Civ/6^ 


Eq .  3 b  reduces  to 
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The  maximum  value  of 


4 


is  determined  by 

* 

maximizing  the  function,  CO'i^ cTTcx. 


elf 


with  respect  to  the  parameter  B  •  A  necessary  condition  for 
the  maximum  is 


Jq  (s/n  coSQ  )  -  0 


?r 


Differentiating  using  the 

product  rule  gives 

4  <St  iA  ^ @  co^>  &  - 

SNA  ^ ©  *  O 

c  os  &  - 

{71 

a  = 

6  3/3  a 

The  fact  that  the  above  value  of  C7  determines  a  maximum 
is  verified  by  the  fact  that  the  second  derivative 

dl  (  SnA  &  coS  o)  -  -  /  O  coSQ  s/  & 
ci&z 


is  always  negative  for  the  admissible  range  of  Q 
Substituting  the  value  of  &  determined  by  Eq.  39  into  Eq.  37 
gives 


io 


wax 
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A- 12  Summary 
1.  The 


£ 

Poo, 


vs.  £  curve  is  constructed  from  the  following: 


CsfCz. 

CA 


■  ®  'j 


t 


]jCt- 


■  i  xk. 

s \*e*  -£■ 


c 

^4 


ifH: !*  cA'^| 


■  '  r 


,  5m  ^  e  §>  ^ 


<-  E 


,  t 


4/ 

41 


where 


v<- ' 


Cl>  -  (\-  O  +  K  (  l  -  § c')  % 

^  -  <  C  1  ‘  O 


<-  i 


xviii 


3.  For  the  linearized  isentropic  case  with  gas  pressures 
in  the  range  15  to  2000  psia,  k  ■  1.40  and  %0  -  .95 
were  used  in  the  calculations. 

4.  For  the  isothermal  case,  k  =  1.0. 


